How precisely will the total cross section be measured at the LHC? 
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It is very likely that hadronic scattering will enter a new regime at the LHC, as the black-disk 
limit is reached. This will lead to a severe change in the t dependence of the real part and of the 
slope of the elastic scattering amplitude, and in turn this may bias the measurement of the total 
cross section. We examine this issue, and suggest new strategies to test the reliability of the total 
cross section measurements. 
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Many models predict that soft interactions will enter a 
new regime at the LHC: given the huge energy, unitari- 
sation may play a crucial role as the central part of the 
protons becomes black. Indeed, all simple-power fits to 
lower energy data will violate unitarity in some partial- 
wave amplitudes before the LHC energy. This means 
that something will happen that will restore it — call 
it saturation, unitarisation or emergence of cuts — and 
this will modify the expectations one has from Rcgge 
models based on simple poles. While conventional mod- 
els predict a total cross section from 90 to 125 mb 
the presence of of a hard pomeron gives around 150 mb 
0, 13 and U- matrix unitarisation can give 230 mb [1] . 
This clearly shows that the uncertainties due to the un- 
derlying models, and especially due to the unitarisation 
scheme, are very large. One may hope to select the true 
models by a measurement of the total cross section. 

The LHC will be well equipped to study in depth the 
diffractive processes, as it will have a superb rapidity cov- 
erage, and two experiments - TOTEM and ATLAS 
Q — plan to measure the total cross section. They intend 
to reach an accuracy on atot of the order of 1%, which 
would indeed give a very stringent test of theory. A few 
assumptions underlie this estimate of the accuracy. First 
of all, the ratio p of the real part to the imaginary part of 
the elastic scattering amplitude, is assumed to be small, 
and to vary little with t: p{s, t) « 0.14 [l|. Secondly, the 
elastic cross section is assumed to fall exponentially with 
t: dN/dt - exp(Bi). 

We want to show here that, if the elastic pp ampli- 
tude reaches a new regime at the LHC, it will invalidate 
the above assumptions and the measurement of the to- 
tal cross section will be biased and much more uncertain 
than foreseen. Indeed, most unitarisation schemes lead 
to novel properties of the elastic amplitude. The prob- 
lems that we talk about here concern a large class of 
models, in which the elastic amplitude contains a fast- 
rising component that needs to be unitarised. We have 
checked that they are present in models saturating the 
profile function, in models using analytic unitarisation 
schemes or in the Dubna Dynamical Model Q . 

To illustrate our point, we shall consider a simple 
unitarised two-component model, which includes a soft 
pomeron and a hard pomeron, and which we shall call 
the eikonalised two-pomeron model (ETPM). This model 



is based on a fit to soft data which includes a hard 
pomeron component @ of intercept 1.4 that accounts for 
the growth of the gluon density at small x in deeply in- 
elastic processes 10]. Although the coefficient of the hard 
pomeron term is small in soft data, it grows like s"'^ so 
that the amplitude will reach the black-disk limit at small 
impact parameter b before the LHC energy 0, [IH . The 
amplitude must then be unitarised and, in this simple 
model, we consider a 1-channel eikonal in the impact- 
parameter representation 0. The net effect of this uni- 
tarisation is to make B{s,t) increase with \t\ at small \t\ 
for LHC energies, as shown in Fig. 1. We also show in 
that figure that the ^-dependence of p{s, t) changes dras- 
tically. We can now explore the consequences of these 
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FIG. 1: Resuhs of the ETPM model: p{s,t) (upper panel) 
and B{s,t) (lower panel) at 100 GeV (plain curve), 500 GeV 
(long dashes), 5 TeV (short dashes) and 14 TeV (dash-dotted 
curve) . 

effects on the experimental measurement itself. We shall 
compare in the following the situation at 2 TeV with that 
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at 14 TeV. We insist that the curves we give for the cross 
sections and for the p parameter are only iUustrative of 
an effect present in many models. The essential ingredi- 
ents are a sizeable value of p, and the strong dependence 
of p and B on t once the black-disk limit is reached. 



I. FITTING PROCEDURE WITH 
LUMINOSITY-DEPENDENT METHOD 

The number of elastic events is related to the total 
hadronic cross section through the following formula: 
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where C is the luminosity, the first term is the Coulomb 
term (a is the electromagnetic coupling constant and 
G{t) the electromagnetic form factor given by G'^{t) = 
{inip - fit) [Arrip - t) / (A^(A - tY) with nip the proton 
mass, A = 0.71 GeV'^ and p, = 2.79), the second term 
is the interference term between the Coulomb amplitude 
and the hadronic amplitude {<t>CN{s,t) is the phase of 
the Coulomb- Nucleon Interference (CNI) term [125), from 
which one can extract p, and the third term is the purely 
hadronic contribution. 

input 





C (fb-1) 


atot (mb) 


p(s,0) 


B(s,0) (GeV-^) 


2 TeV 


1 


82.7 


0.23 


18.7 


14 TeV 


1 


152.5 


0.24 


21.4 


output for C and p fixed 


2 TeV 


1 


83.61 ± .44 


0.15 


23.6 ±0.2 


2 TeV 


0.95 


85.8 ± .45 


0.15 


23.6 ±0.2 


output for all parameters free 


2 TeV 


0.93 ± 0.07 


85.2 ±3 


0.15 ±0.04 


18.1 ±0.25 


14 TeV 


1.15 ±0.05 


142.3 ±2.8 


0.29 ±0.06 


23.6 ±0.2 



TABLE I: Input parameters for the simulated data at y/s = 
2 TeV and = 14 TeV obtained in EFTM model, and the 
results of fits to these data with a simple exponential form of 
the scattering amplitude. 
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FIG. 2: The simulated data at = 14 TeV, the theoretical 
curve from which the data are generated (plain) and the fit 
to them for p fixed at 0.1 (dashed). 
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FIG. 3: The size of atot obtained by fitting the simulated data 
assuming a fixed value of p at ^/s — 2 TeV with £. fixed at 
the input value (upper panel) and with free C (lower panel); 
the diamond gives the true value of atot and p{s, 0). 



We can use this formula to generate simulated data for 
two energies -y/s — 2 TeV and y/s — 14 TeV, using B{s, t) 
and p{s,t) calculated in the ETPM. We assume that 90 
points will be measured in a i interval identical to that of 
the UA4/2 experiment -0.1 GeV^ < t < -0.0006 GeV^ 
[l^. We then randomize the theoretical ("true") curve 
assuming Gaussian errors similar to those of UA4/2. The 
resulting simulated data are shown in Fig. 2 for y/s = 14 
TeV, and correspond to the parameters given in Table 1. 



One can then fit these simulated data according to Eq. [T] 
but assuming constant B and p. One has 2 extra param- 
eters besides p and B: £, the luminosity, and the total 
cross section atoti which is what one aims to measure. 

The result of the fitting procedure at 2 TeV is shown 
in Fig. 3, where the correlation between the value of p 
and atot is shown. If C is fixed at 1 fb~^ {i.e. if we know 
the luminosity), the difference between the central values 
of the fitted atot is small « 0.3 mb; the errors from the 
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fitting procedure are 1 mb and the obtained value of atot 
differs from the input by 0.5 mb. A very different picture 
appears if we fit the luminosity. We obtain then an error 
^(^tot — 1-9 mb and the correlation between atot and the 
assumed value of p is very high, as seen in Fig. 3 and 
Table 1. The result of a joint fit to C, B, atot and p is 
also shown in Table 1. All parameters are within 1 a, 
except p(0) which is 2a from the input. 

The situation at ^/s — 14 TeV is much worse: one sees 
from Fig. 4 that the result of the correlation between 
p and atot increases drastically: even if one knows the 
luminosity, the dependence of atot on p is very strong. 
The difference between the true atot and the fitted value 
reaches 2.5 mb, as shown in Table 1. All the parame- 
ters are now several standard deviations from their true 
value. Note also that, because the CNI term is negative, 
a decrease in p from its true value 0.24 to 0.1 leads to an 
increase in the value of da/dt, as seen in Fig. 2. 



if one allows p to range from 0.05 fl*] to 0.3, one will get 
a 10 % change in atot, which will only be added to the 
uncertainty coming from the estimate of N^i + Ninei ■ So 
it seems to us that it is illusory to hope for an accuracy 
on atot better than 10 % from this method. 
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II. FITTING PROCEDURE WITH THE 
LUMINOSITY-INDEPENDENT METHOD 

Another way to extract the total cross section (for ex- 
ample, see [iJl) is the luminosity- independent method, 
which gives 



167r {dNei/dt)\t=o 



(2) 



To obtain the information needed for this method one 
must measure the elastic rate at values of |i| large enough 
to neglect the Coulomb amplitude Hence, one hopes 
to obtain a more accurate result with less information. 

This method relies on the hope that N^i + Ni^ei can be 
measured accurately. However, there are three problems. 
The first one concerns the Coulomb and CNI regions: one 
needs to cut them off, but Nei + Ninei is the total number 
of hadronic events, so one must compensate somehow for 
that cut. The second problem comes from the fact that 
part of the inelastic events (such as N* production) will 
escape the detector. The third problem, which we shall 
address here, concerns the extrapolation of dN/ dt from 
a minimum \t\ far away from the CNI region to i = 0. 

To simulate this analysis, we take iVg/ -I- Nin — natot, 
and let n go from 0.9 to 1.1. We adopt the cuts planned 
for the TOTEM experiment at the LHC, -t in [0.03, 0.1] 
GeV^ , so we are left with in 51 simulated data points. 

For the analysis at 2 TeV, if we take n = 1 the depen- 
dence of atot over p will be comparable to that shown 
in Fig. 2, but the errors on atot will increase. For ex- 
ample, if we fix p = 0.15, then atot — 83.61 ± 0.44 mb. 
atot changes by 1.9 mb when p goes from 0.05 to 0.25. 
If n can go as low as 0.95, then the measurement of atot 
increases, e.g. atot — 85.8 mb when p = 0.15. 

At i/i = 14 TeV, the changes are more pronounced, 
as seen in Fig. 4. Cutting off the CNI region removes 
the possibility to measure p. As the normalization is 
inversely proportional to 1 + p'^, it is rather obvious that 



FIG. 4: atot at — 14 TeV with fixed normalization (n — 1) 
(central curve) andn = 1±0.1 (exterior curves). The diamond 
indicates the input atot and p(s, 0) from which the data were 
simulated. 

The large discrepancy between the measurements of 
atot by CDF [l| and E710/E811 [3 probably has its 
origin in the uncertainty on A^e; + A^in and on the de- 
pendence of p and B on t, and not in some experimental 
mistake. Hence this large difference really reflects the 
real error one is to expect from this method. 



III. CONCLUSION 

Our analysis shows that both methods introduce large 
correlations between p and atot- As it is very likely that 
at the LHC unitarisation will play an important role, 
one should not assume that B[s,t) and p{s,t) are con- 
stant with t, and their exact behavior with t is model- 
dependent. Hence the inescapable conclusion is that a 
1 % measurement of atot will only be possible if one mea- 
sure B{s,t) and p{s,t) as well. 

\t\ should go from very small values, as close as pos- 
sible to zero, to about 0.1 GeV^, with sufficiently small 
bins. And it will be important to allow all parameters 
to vary. Indeed, the measurement of p performed by 
UA4/2 (0.135 ± 0.015) [13] seemed to contradict that of 
UA4 (0.24 ± 0.02) ^ only because atot was fixed. Al- 
lowing atot to be fitted to leads to an agreement between 
the two measurements [l7| . 

We also believe that the luminosity-dependent method 
is preferable, as it uses more information. The mea- 
surement of p performed at the Tevatron [3] used the 
luminosity-independent method with very large bins in 
t, and the interval considered was 0.00095 < \t\ < 
0.1431 GeV^. On the lower side, one reached very small 
t, so that the behavior of the amplitude cannot be taken 
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as a single exponential because of the CNI effect. On the 
upper side, the intervals in t were too big to measure the 
specific properties of the CNI region. Hence the p pa- 
rameter extracted is very uncertain, and it could be that 
it varied appreciably with t. 

As the standard fitting procedure can give misleading 
results, we remind the reader that additional methods 
have been proposed to check the validity of the assump- 
tions entering the fits. Firstly, it is possible to extract 
the value of p{s,t) at small but nonzero t [H, HO]) using 
the fact that the Coulomb amplitude Fc{t) has an op- 
posite sign to that of the real part of the pp amplitude 
Fn{t), so that there is a value t = tc for which Fc{tc) — 
—^eFn{tc), so that the differential cross section da/dt 
has a local minimum at that value. The position of 
this minimum depends strongly on the form assumed for 
p{s, t) and extracting its value would show whether p 
varies quickly with t or not. Using this method, it was 
found in [19, 20] that already at ^/s — 52.8 GeV, p was 
not constant with t. 

It is also possible, at small to determine the elas- 
tic hadronic cross section via an iterative method, which 
takes advantage of the expression of the total elastic cross 
section in the CNI region [2l| . 

Finally, one can adapt a method that was first designed 
to study eventual oscillations in da/dt ^2^. The idea is 
to compare two statistically independent samples built 



by binning the whole t range in small intervals, and by 
keeping e.g. one interval out of two. The deviations 
of the experimental values from theoretical expectations, 
weighted by the experimental error, are then summed for 
each sample k: 

AR\t)= ^ [(daVrfi.r^ - (d^M.)"']/<5^^ (3) 

|td<l*l 

where S^^^ is the experimental error. If the theoretical 
curve does not precisely describe the experimental data, 
(for example, if the physical hadron amplitude does not 
have an exactly exponential behavior with momentum 
transfer), the sum Ai?'^(i) will differ from zero, going 
beyond the size of the statistical error. 

Using these methods will help to test the assumptions 
entering the future experimental analyses of TOTEM and 
ATLAS, and may lead to a much more reliable measure- 
ment of p and atot- 
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